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ABSTRACT 

During the observation period 1997, the nearby Blazar Mkn 501 showed extremely strong 
emission and high variability. We examine multiwavelength aspects of this event using radio, 
optical, soft and hard X-ray and TeV data. We concentrate on the medium-timescale variability 
of the broadband spectra, averaged over weekly intervals. 

We confirm the previously found correlation between soft and hard X-ray emission and the 
emission at TeV energies, while the source shows only minor variability at radio and optical 
wavelengths. The non-linear correlation between hard X-ray and TeV fluxes is consistent with a 
simple analytic estimate based on an SSC model in which Klein-Nishina effects are important 
for the highest-energy electrons in the jet, and flux variations are caused by variations of the 
electron density and/or the spectral index of the electron injection spectrum. 

The time-averaged spectra are fitted with a Synchrotron Self-Compton (SSC) dominated 
leptonic jet model, using the full Klein-Nishina cross section and following the self-consistent 
evolution of relativistic particles along the jet, accounting for 77 absorption and pair production 
within the source as well as due to the intergalactic infrared background radiation. The 
contribution from external inverse-Compton scattering is tightly constrained by the low 
maximum EGRET flux and found to be negligible at TeV energies. We find that high levels of 
the X-ray and TeV fiuxes can be explained by a hardening of the energy spectra of electrons 
injected at the base of the jet, in remarkable contrast to the trend found for 7-ray fiares of the 
flat-spectrum radio quasar PKS 0528-1-134. 

Subject headings: BL Lacertae objects: individual: Mkn 501 — radiation mechanisms; 
non-thermal 
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Introduction 



The BL Lac object Mkn 501 is very close {z = 0.0337, Ulrich et al. 1975) and has been studied 
extensively at all wavelengths. Together with its sister object Mkn 421, it was among the two first BL Lac 



objects with known radio (CoUa et al. 1972), X-ray (Schwartz et al. 1978) and TeV gamma-ray (Quinn et 



al. 199£, Bradbury et al. 1997) counterparts. Recently, it was also marginally detected at photon energies 
> 100 MeV by the EGRET instrument on board CGRO (Kataoka et al. |l999| ). 



During 1997, the object was found to be in an extreme high state with a TeV flux on average 20 times 



higher than in 1996 (Breslin et al. 1997). The source exhibited strong variability on timescales of days with 
a possible quasi-periodically varying component with a timescale of about 25 days (Kranich et al. 



19991) 



To complete the list of reasons for excitement for the observers, Mkn 501 reached in some of its flares fluxes 



of more than 10^"'^'^ cm~^s~^ (above 1.5 TeV) - the most intense TeV emission ever measured so far from 
any astronomical object. However, the shortest observed variability timescale (5 hours, Aharonian et al. 
1999a) was significantly longer than that observed for Mkn 421 (Gaidos et al. 1996| ). 



1998, 



All available TeV observatories monitored the event for several months (Samuelson et al. 

1999| ). The most complete 



Aharonian et al. 1999a 1999b, Hayashida et al. 1998, Djannati-Ata'i et al 



dataset was produced by the HEGRA Cherenkov Telescope "1" (CTl) (Aharonian et al., 1999bD which was 
even able to observe Mkn 501 under the presence of moonlight, however with reduced sensitivity, thereby 
filling many gaps in the lightcurve. This telescope also obtained the confirmatory observations in 1996. 

Also from HEGRA come probably the most accurate spectral measurements in the TeV regime. They 
were carried out by the HEGRA system of (at the time 4) Cherenkov telescopes (GTS) and are largely 



concurrent with the CTl measurements, however with less time coverage (Aharonian et al. 1999a 



The origin of the TeV 7-ray emission and the reasons for its variability are still essentially unknown. 
The most popular models explain the TeV emission as near-infrared to UV photons which have been 
upscattered via the inverse Compton effect by very high energy electrons which are known from radio, 
optical and X-ray observations to be present in the jets of BL Lac objects. Possible sources of the 
seed photons for Compton scattering are the synchrotron radiation produced within the jet by the 
same population of relativistic electrons (synchrotron self-Compton model, SSC; Marscher & Gear 1985| ; 



Maraschi, Celotti & Ghisellini 1992; Bloom & Marscher 1996), or radiation from outside the jet (external 
inverse Compton model, EIC). This external radiation could be the quasi-thermal radiation field of an 
accretion disk surrounding a supermassive black hole which is generally believed to power the relativistic 
jets. The accretion disk radiation can enter the jet either directly (Dermer, Schlickeiser & Mastichiadis 
1992| , D ermer & Schlickeiser |l993| ) or after being rescattered by circumnuclear material (Sikora, Begelman 



& Rees 1994 , Blandford & Levinson 1995 , Dermer, Sturner & Schlickeiser 1997 ). It is also possible that 
synchrotron radiation produced within the jet and refiected by circumnuclear debris is the dominant source 
of soft photons during fiares (Ghisellini & Madau 1996, Bednarek |1998| ), although it has been shown that 
this process is unlikely to be efficient in the case of BL Lac objects (Bottcher & Dermer 1998). 



As an alternative, Mannheim ( 1993 ) has suggested a hadronic model in which protons are the primarily 
accelerated particles in the jet and the 7-ray emission is produced by secondary pions and electron-positron 
pairs produced in photopion and photopair production interactions of the ultrarelativistic protons in the 
jet with external radiation. The time-averaged broadband spectrum of Mkn 501 has been fitted using this 



model by Mannheim et al. (1996, 1998). However, the attempts to explain the short and intermediate-term 
variability of blazars with hadronic jet models have only just started (Rachen & Mannheim 1999| ). For this 
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reason, we concentrate on leptonic jet models in this paper since the variabihty time scales predicted by 
these models are in good agreement with the observed intraday variability of blazars. 



Recently, Fossati et al. (1997) have compared the broadband spectra of different types of 7-ray emitting 
AGN and suggested a continuous sequence FSRQs (flat spectrum radio quasars) LBLs (low-frequency 
peaked BL Lac objects) HBLs (high-frequency peaked BL Lac objects), characterized by decreasing 
bolometric luminosity, decreasing dominance of the total energy output in 7-rays compared to the emission 
at lower frequencies and a shift of the peak frequencies of the synchrotron and the 7-ray component towards 
higher frequencies. Recent modelling efforts of various blazar-type AGN have revealed that this sequence 
is consistent with a decreasing importance of external radiation as a source of soft photons for Compton 



scattering in the jet (Ghisellini et al. 1998| ). This suggests that the extreme HBLs like Mkn 501 or Mkn 421 



can be well fitted with strongly SSC-dominated jet models, as was shown, e. g., by Mastichiadis & Kirk 



(1997) and Plan et al. (199S 



From the SSC model one expects a strong correlation between the X-ray and the TeV emission. And 
indeed, by comparing the daily TeV measurements with daily averages from the All Sky Monitor (ASM) on 
board the Rossi X-Ray Timing Explorer (RXTE) , a significant correlation with a most probable time-lag of 



days (no time-lag) was found (Aharonian et al. 1999a , 1999b ). Furthermore, a clear overall X-ray high 



state was visible in the 1997 ASM measurements which coincided with the TeV high state (see also figure 
|l|). And, as discovered by BeppoSAX, the X-ray peak of the spectral energy distibution (SED) had shifted 
from 10-20 keV in 1996 to 100-200 keV in 1997 (Pian et al. |1998| ), consistent with the assumption that the 



X-ray and TeV-7-ray flares are produced by a more powerful electron acceleration, shifting the high-energy 
cutoff of the electron distribution to higher energies and hardening the X-ray spectrum. This observation 
was, however, only made during one flare in 1997. In this paper we give evidence that the hardening took 
place during all flares in 1997. An alternative explanation for the synchrotron spectral changes, in terms of 
a steadily-emitting helical-jet model, has been presented by Villata & Raiteri 



Apart from the short-term variability on timescales of hours to days, there is also a longer-term 
variability in Blazars which has so far been investigated mainly in the optical regime (see e.g. Katajainen 
et al. ( 1999 ) and references therein) and recently also in X-rays (e.g. Mc Hardy 1999| ). This variability 



shows remarkable amplitudes (e.g. 4.7 mag for Mkn 421 in the optical) and in the case of OJ 287 there is 
even evidence for periodicity. 

In this respect Mkn 501 is not yet very well explored. Mkn 501 has been observed in the TeV regime 
since its discovery in 1995. The source showed low emission close to the sensitivity limits until the onset 
of the 1997 high state. Since the duty cycle of the TeV observatories is only about 10 %, this docs not 
prove the absence of strong short flares prior to this high state, but the increased average intensity of the 
source can be described as an increase in flaring probability from 1996 to 1997 by at least an order of 
magnitude. This description is especially appropriate since even during the high state, the source returns 
to quiescent (comparable with 1996) levels of emission for periods of up to a few days. These transitions 
are seen irrespective of the presence of short, strong flares on time-scales of a day. 

The last observations in 1996 were made in August and found Mkn 501 still quiescent while the flrst 
observations in 1997 were made in February and found the source already flaring. The transition from low 
to high flaring probability obviously took place on timescales shorter than half a year. The correlated X-ray 
data from the RXTE ASM confirm this and show in addition that the change has been a smooth process 



over several months (see e.g. figure 8 in Aharonian et al. 1999b) 



In this article we explore the multi-wavelength variability of Mkn 501 over medium timescales, and try 
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to relate its behaviour to the physical parameters of a leptonic jet model. For this purpose we construct 
weekly SEDs using the HEGRA CTl flux data and HEGRA CT System spectral data together with data 
from longer wavelengths, namely radio, optical, soft X-ray and hard X-rays. See table ^for an enumeration 
of the instruments and energy ranges. 

This paper may represent an important step toward understanding the origin of strong changes in the 
flaring probability of Mkn 501 and Blazars in general, but our data is clearly not sufficient to give a final 
answer to this question. 

All data except that from HEGRA and the RXTE ASM are published here for the first time. The 
BATSE data are especially valuable since they confine the intensity at the X-ray peak of the SED. 



2. Observations and Data Analysis 

The observations used to fit the multi-wavelength spectrum cover mostly the synchrotron part 
of the SED. Only the TeV data explore what is believed to be the Inverse Compton emission, 
although OSSE and EGRET observations cover a few days in 1996 and 1997. As a guide line for 
our fitting procedure we take into account the highest ever observed EGRET flux from Mkn 501, 
F{> 100 MeV) — (32 ± 13) • 10^^ photons cm^^ s^^ with a photon spectral index a — 1.6 ± 0.5, measured 



during 1996 Mar 25 - 28 (Kataoka et al. 1999), as an upper limit 



In the definition of time-bins for the multi-wavelength dataset, we start by subdividing the HEGRA 
CTl lightcurve and bin all other data accordingly. 

The data from HEGRA are binned in time from March through October 1997 (MJD 50514 - MJD 
50708) resulting in 28 equidistant time bins. This weekly temporal resolution is an order of magnitude 



larger that the longest observed TeV intraday variability timescale of 15 h (Aharonian et al. 1999a) and 
is believed appropriate, given the nature of the available data and the medium scale variability timescale 
which we have chosen to explore. 

However, the HEGRA points are not spread as uniformly over time as are the BATSE and RXTE 
data. The "center of gravity in observation time" (defined as the weighted mean of the observation time of 
the daily points each weighted by the duration of the individual observation) from HEGRA would therefore 
in general not coincide with that from BATSE and RXTE. In order to compensate for this, we calculate the 
center of gravity in observation-time for each weekly HEGRA timebin and use these as time-bin- centers for 
the other data. The edges of these bins are then defined by the average of two adjacent bin-centers. 

Figure |l] shows the lightcurves from each instrument which went into this analysis. The weekly points 
after the binning are given in table |[ The first column in the table gives the bin-centers as described above. 
The following subsections describe the data used to compile the table. 



2.1. Optical and Radio Observations 



The optical observations were performed using the telescopes and filters listed in table |^. All the 
observations were made with CCD-cameras. All CCD-images were treated the normal way with flat field 
and bias corrections. The magnitudes were measured using either DAOPHOT in IRAF (NOT and Tuorla 
data; for more details see Katajainen et al. 



199£), or with the ROBIN-procedure developed in Torino 
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(e.g. Villata et al. 1997), or using the automated reduction routine developed in Perugia (e.g. Tosti et 
1996). All magnitudes were measured using 10 arcsecond aperture. The use of the same aperture is 



al 



important because Mkn 501 has a large host galaxy (e.g. Nilsson et al. 1999 ). For the calibration we used 
the calibration sequence given by Fiorucci & Tosti ( |l996| ) and Villata et al. ( |l998D . 



For some of the weekly observing periods, no optical data were available. However, the optical flux is 
extremely important for constraining the spectral index of the synchrotron spectra in the optical - X-ray 
range, which, in turn, is essential for constraining the spectral index of the electron spectrum. As can 
be seen in table |^, the optical flux exhibits only moderate variability on the time scales considered here. 
Thus, for fitting purposes we assume that the optical flux during those viewing periods for which no optical 
observations were available, was within the range of optical fluxes observed during the whole campaign 
which yields uF^{opt. 



(7.485 ± 0.774) • 10"" erg cm'^ s" 



The radio observations were made with the Metsahovi Radio Telescope at 22 GHz as part of the 
ongoing quasar monitoring program, which started in 1980. Currently about 85 sources, mainly Northern 
flat spectrum quasars, should be observed monthly at 22 and 37 GHz. The total number of observations is 
now over 40000. For details of the observing strategy and reductions, as well as the data until 1995.5 see 
Terasranta et al. ( |l998| ). 



The integration time of each point was essentially the same, so single points were formed from the data 
points in the individual timebins by simple averaging. 



2.2. Soft X-Ray Data 



Since the beginning of 1996, the all-sky monitor (ASM) of the Rossi X-ray Timing Explorer (RXTE) 
satellite has been observing Mkn 501 in the 2-12 keV energy band. The data used in this analysis are taken 
from the publicly available ASM data products provided by the ASM/RXTE teams at MIT and at the 
RXTE SOF and GOF at NASA's GSFC. These measurements are given by the authors (see e.g. Levine et 



al. 1996) as rates i?i,i?2,i?3 in units of counts per second in three energy bins: 1.3-3.0 keV, 3.0-5.0 keV, 
and 5.0-12.1 keV. From these we calculate the total rate R by summing the three bins. 



For Mkn 501 we obtain by averaging over the period MJD 50510 - 50710 



R 



1.3 ±0.4 



(1) 



In order to assess the X-ray spectrum and flux of Mkn 501, we use the ASM data for the Crab Nebula 
which is publicly available from the same source. The average rate of the Crab Nebula during the period 
MJD 50510 - 50710 is 

i?c™6 = 75.7 ±0.5 (2) 

Hence, in soft X-rays, Mkn 501 is even during this high state a significantly weaker source than the Crab 
Nebula. 

In the energy range 2 keV to 60 keV, the spectrum of the Crab Nebula is a stable power law. The 
spectral index of the differential photon spectrum is known with an accuracy of 1% to be acrab = 2.1 (e.g. 



Toor & Seward |1974| , Pravdo & Serlemitsos |1981| or Pravdo, Angelini & Harding |1997| ). The emission 
measured by the RXTE ASM is the sum of the steady Nebula and the pulsed Crab Pulsar emission. The 
latter has on average a harder spectrum than the nebula. Below 12.1 keV, however, the pulsed fraction of 
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this emission is < 8 %. We can therefore use the knowledge of the Crab Nebula spectrum to derive an 
approximate relative calibration factor k for the rates R2 and i?3 using 

— = ^ ^ = 1.213 (3 

R3,crab i;!-"-"' - V ; 

where E2 = 3.0 keV, = 5.0 keV, and E4 — 12.1 keV are the energy bin edges of the second and third 
energy bin. We ignore the first energy bin since it is strongly influenced by interstellar X-ray absorption 
which is dependent on the column density and hence varies between sources (Remillard |l999D . 



Averaging over the available Crab data from the observation period under discussion (MJD 50510 - 
50710) we obtain k — 1.318 ± 0.0024. There is no indication of a variability of the value of k (see figure |]). 
Hence we can assume that it is also valid for the Mkn 501 observations of the same detector. In order take 
into account that there is a pulsed component with a harder spectrum, we add an additional error of 8 % 
to the error of R3 and obtain thus 

yfc = 1.32 ±0.023 (4) 



For Mkn 501 we calculate the spectral index a of the differential photon spectrum using the equation 

?2,mfcn501 _ E^ 



kR2,mkn501 _ -E3 " — £^2 " 
R3,mkn501 E\^°' — 



(5) 



and varying a until the two sides of the equation are equal to an accuracy better than 0.1 %. The error 5a 
of this index, we estimate from the approximate formula 

^ log(fc ■ i?2/(5.0 - 3.0)) - log(fl3/(12.1 - 5.0)) 
" log((3.0 + 5.0)7(5.0 +12.1)) ^ ' 

which leads to 



fo= ^1.733. ((|:)^ + (^)^ + (^)^) (7) 

where 5fc, (5i?2, and SR^ are the errors of the corresponding quantities. The time resolved values for amknboi 
are shown in figure |. The same method applied to the Crab Nebula data (however with 6k = 0.0024) 
yields, as expected, on average the spectral index we have put in. Also this is shown in figure ^. From a 
constant fit to the values in this figure, we obtain 

amfc„50i = 1.76 ±0.04 (8) 

and 

ae™6 = 2.10 ±0.02 (9) 

Thus we find that during the 1997 high state, Mkn 501 had on average a significantly harder soft X-ray 
spectrum than the Crab Nebula. The spectral variability of Mkn 501 is below the ASM's sensitivity. The 
distribution of the weekly values is consistent with a constant value (reduced = 0.93) and so is that for 
the Crab Nebula (reduced = 0.95). 

In order to calculate the vF^ values in erg cm~^ s""'^, we use the knowledge of the fiux of the Crab 
Nebula. Here we face the problem that there is still a disagreement of up to 25 % between the measured 
normalization constants of the Crab Spectrum from different experiments although there is perfect 
agreement in the spectral index. The discrepancy seems to stem from not very well understood systematic 



differences between the detectors (see the discussion in Pravdo, Angelini & Harding 1997). We use the 
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two most extreme recent measurements of the differential photon flux of the Crab Nebula together with 
the Crab Pulsar (pulse-averaged) and take their average as our normalization and their difference as the 
systematic error of this quantity. From Pravdo & Serlemitsos ( 198f ) we get at 5.2 keV a flux of (0.236 
± 0.006) photons cm-2s-ikeV"\ while from Pravdo, Angelini & Harding ( |f997[ ) we get (0.302 ± 0.001) 
photons cm^^s^^keV^^. The average of these values corresponds to a differential energy flux of 

i^crafc (5.2keV) = (2.24 ± 0.27syst) x lO'^erg cm-^s^^keV^^ (10) 

where the systematic error is the difference between the averaged values divided by 2. The statistical errors 
of the two measurements are negligible. The energy flux of any other source with ASM rates i?2 and R3 
and differential photon spectral index a is then calculated by 

„p lay^nr^-'^^-i] - (R mRW^-^] Fgrab (5 . 2keV) [crg cm~^s~^keV~^] 
i/i'i.[ergcm s J - (H2 + n3)[s \- — rr-^ ^ 

(,-n,2,crab + tv3,crab j [S J 

.5 2i+"c™b-a . 1-a . £4 (11) 

where we correct for the difference in the spectra of Crab Nebula and the source in question by making the 
Ansatz that the measured differential rates have the same ratio as the differential fluxes. Furthermore we 
have used the result i?2,cra6 + Rz,crab — 48.5 ± 0.7 s^^ obtained from the dataset under discussion. Inserting 
the spectral index of Mkn 501 gives: 

z/i^^(Mkn501,5.2keV) = (R2 +R3)[s-^] • 2.40 x 10-^°erg cm'^g-i (12) 

The error of this flux value is determined by propagating all errors of the quantities involved which gives 

(5(iyF^)(Mkn501,5.2keV) 



(v/5.76 X 10-20 . {{5R2Y + {SRzf) + 1.21 x 10-23 . {R^ + i?3)2 + 0.12 • i/i^^) erg cm-2s-i 
where the term outside the square-root stems from the systematic error of the Crab flux normalization 



(13) 



2.3. Hard X-Ray Data 

The Hard X-ray fluxes were measured using the Burst and Transient Source Experiment (BATSE) 
onboard the CGRO satellite. Although BATSE is an uncoUimated detector, accurate point source fluxes 



can be measured using the Earth Occultation method described in Harmon et al. (1992). Daily sensitivities 
are 100 mCrab and over integrations of years, sources as weak as 3 mCrab can be detected. Several months 
are usually needed to obtain a statistically significant flux from Mkn 501 with BATSE, but intense daily 
flares can be seen, and the overall high state of 1997 allowed useful measurements in each weekly interval 
even outside the flares. 

The source is measured only when it sets and rises from behind the Earth. Since two such occultations 
occur per spacecraft orbit (roughly every 90 minutes), up to 32 independent flux measurements can be 
made per day. There exists, however, a wide variation in this number because of passage of the spacecraft 
through the South Atlantic Anomaly, telemetry gaps from loss of TDRSS contact, and other events, which 
occur randomly relative to the Mkn 501 steps. For the data shown here between 61 and 211 measurements 
went into an individual weekly point. Each measurement lasts about 8 seconds giving a duty cycle of about 
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0.2%. In a 7 day time bin during the period discussed in this paper, BATSE observes of the order of 10^ 
photons from the source (more during the flares). 

The fluxes are integrated between 20 and 200 keV, and are calculated by folding the measured counts 
through the BATSE detector response assuming a differential source powerlaw spectrum of index -2.0. The 
-2.0 spectrum is the best fit to the flare measured on MJD 50550-51 (between 20 and 1000 keV). Spectra for 
other time intervals were also calculated and were consistent with -2.0. Uncertainties of 10 % in the spectral 
index during flare times, larger outside, make spectral variability difficult to assess for this source, so that 
the index of -2.0 was used for each weekly interval. Fluxes were calculated for smaller energy bands but the 
single 20-200 keV (median energy 36.4 keV) point for each interval is the most useful outside intense fiares. 

Seven of the weekly averages are not statistically significant, and one period shows a 1.8ct deficit. These 
eight points arc inconsistent with a zero-level flux at the 99 % confidence level. The seven low but positive 
points are inconsistent with a zero-fiux level at the 98-99 % confidence level and are well fit by a constant 
flux of 1.1 ± 0.3 X 10^^" erg cm~^ s~^, implying that a flux is present which is below the sensitivity of the 
BATSE instrument in weekly integrations. An analysis of 30 blank fields on the sky shows that with 100 
weekly integrations for each field one can not distinguish (using an F-test) between a zero-flux level and the 
weighted mean of the weekly averages for any of the 30 fields. This implies that systematic effects are not 
responsible for the excess seen in the low points and that a flux is indeed present. For this reason these low 
but positive values have been included as detections in the multiwavelength fits. 



2.4. TeV data 

As described above, we use the data from the CTl lightcurve of Mkn 501 (integral flux above 1.5 TeV) 



published by Aharonian et al. ( 1999b ). These data are available in daily points based on observations 
of between 0.5 h and 5 h duration each. We group these points according to our weekly time-bins and 
calculate an average flux from the up to seven values weighting each daily point by its observation time. 



Apart from giving daily flux measurements, the HEGRA papers Aharonian et al. ( 1999a ), (1999b) 



and (1999c) also determine the average spectral shape in the range between 0.5 and w 25 TeV with high 
accuracy. They find 

dF/dE = NoE-°' cx-p{-E/Eo) (14) 

where No = (10.8±0.2±2.1) x 10"" cm^^ TeV^^, a = 1.92±0.03±0.20, and Eq = (6.2±0.4±2.2) TeV. 
The first error given is the statistical, the second the systematic error. Furthermore, they find that there is 
no spectral variability up to their sensitivity of 5a < 0.1 on all relevant timescales. 

In order to include this important spectral information in our model fit, we make the assumption that 
there is indeed no spectral variability and extrapolate the points measured at 1.5 TeV by CTl up to 10 TeV 
and down to 0.8 TeV. From the average integral flux values in units of photons cm^'^s^^ we obtain the 
i^F^ values at photon energy E in TeV using 

.FAE)[eTg cm-^s-i] = 1.6022 • F,., g'""^ ' '^^^'^/^^^ . (15) 

/ dee-" exp(-e/Eo) 
1.5 

The extrapolation uses the measured spectral shape (equation |l^) and fully propagates all statistical errors 
(error of the CTl point, error of the spectral index a and the error of the cut-off energy Eq) to form the 
error of the extrapolated point. The systematic errors of the spectral shape are expected to influence all 
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measured TeV spectra in the same way since they are caused by the 15 % uncertainty in the absolute energy 
cahbration of the Cherenkov telescopes. Thus, for the purpose of a spectral variability study, the statistical 
errors are those which actually determine the uncertainty in the differences of the spectral shape between 
different measurements. For this reason we propagate only the statistical errors in the extrapolation. In 



this way, the columns 6 and 8 of table are obtained. In the model fit (section 2.5) we take into account 



the uncertainty in the energy scale by introducing errors of ±15 % along the energy axis. 

The energies to which we extrapolate are chosen as a compromise of maximum distance to 1.5 TeV and 
minimum systematic errors. The latter are increasing up to several 10 % towards both ends of the range for 



which the spectrum has been measured (see Aharonian et al. 1999c, figure 9) but are still small at 0.8 and 
10 TeV. 

The correlation between the three TeV points which we introduce into the model fit by performing the 
described extrapolation is not problematic since we do not plan to calculate absolute values in the fits 
for proving that the model describes the data better than another. Instead, the fits serve the aim to study 
the time-dependent behaviour of the model parameters. The extrapolated points are only a means to take 
into account the available spectral information. 



2.5. Fitting the model to the SEDs 



To each weekly SED we fit the Blazar jet model described in detail by Bottcher et al. ( 1997] ). The model 



assumes that isolated components (blobs) of relativistic pair plasma, which are assumed to be spherical 
in the co-moving frame, are injected instantaneously into the jet, and follows the self-consistent evolution 
of the particle and radiation spectra as the blob moves outward along the jet, taking into account all 
relevant radiation, cooling, and absorption mechanisms: synchrotron radiation, synchrotron self-absorption, 
synchrotron self-Compton scattering, external Compton scattering of direct accretion disk radiation, 7-7 
absorption and pair production intrinsic to the source. The magnetic field is assumed to decay along the 
jet as i? oc r~^, where r is the distance from the center of the AGN. The emerging, time-averaged spectra 
are corrected for 77 absorption by the intergalactic infrared background radiation using the lower model 



spectrum given by Malkan & Stecker (1998) 



Since we are interested in weekly averages, the emission from single blobs is time-integrated over the 
jet evolution and subsequently re-converted into a flux by dividing the fluence by an average repetition 
time Atiop of blob ejection events. We assume that a fraction / < 1 of the jet is filled with relativistic pair 
plasma. The filling factor is given by / '--^ R'g / (T c At,-cp) , where R'g is the blob radius in the co-moving 
frame and T is the bulk Lorentz factor of the blobs. Even if the filling factor is close to unity, it is a 
reasonable approximation to assume that the blobs do not interact with each other because the synchrotron 
and SSC radiation produced within the jet are isotropic in the comoving frame so that most of the radiation 
escapes to the sides without interaction with the rest of the jet. 

As mentioned in the introduction, extreme HBLs like Mkn 501 or Mkn 421 are generally well described 
by a pure SSC model. A simple analytic estimate shows that the observed TeV 7-ray spectrum can 
not plausibly be produced by Comptonization of radiation from an accretion disk around the putative 
supermassive black hole in the center of Mrk 501: The spectrum emitted by an optically thick, geometrically 
thin accretion disk is reasonably well approximated by a blackbody spectrum whose temperature, for a black 
hole mass of > 10® Mq, yields an average photon energy of the disk radiation of ejj = hvo / {m^c?) ^ 10^^. 
If external Comptonization is to be efficient in competition with the synchrotron sclf-Compton process, the 
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blob has to be rather close to the accretion disk, or a significant fraction of the disk photons has to be 
rescattered into the jet by surrounding material, so that the bulk of disk photons enters the blob from the 
side and is blue shifted by a factor of ^ F into the blob rest frame. Thus, due to the strong reduction of the 
Klein- Nishina cross section for 7ee' > 1 (where e' is the photon energy in the comoving frame), no significant 
radiative output at (observer's frame) energies Cots > D/{eD^) ^ 10^ (where D is the Doppler factor), 
corresponding to Eots ^ 100 GeV, will be produced by external Comptonization, ruling out this process 
to explain the observed high-energy spectrum extending to TeV energies. If the high-energy spectrum in 
the ^ 100 GeV regime is produced by the SSC process, then the level of SSC radiation at 1 GeV may be 
estimated by 

p-3 

.F^^^^il GeV) « .F^'^^iE,,) ' (16) 

where Epk ^0.1 TeV is the energy of the vF^ peak in the high-energy part of the spectrum and p is the 
spectral index of the injected electron spectrum. Inserting a typical value of vFv{Epk) ^ 10^^ JyHz and 
p ~ 2.5, this yields vF^^'-^ [1 GeV) ^ 3 • 10^^ JyHz, which is already close to the maximum ever observed 
of ~ 10^^ JyHz. In reality, the SSC spectrum is not a straight power-law below Epk, but shows a gradual 
turnover so that the actual SSC flux at 1 GeV is substantially higher than the above estimate, leaving little 
room for an additional EIC contribution. For our fitting procedure, we are thus using a simplified version of 
our jet simulation code, in which the photon output (but not the electron cooling) from Comptonization of 
accretion disk radiation is neglected. These simulations properly account for the self-consistent cooling of 
the electron population and 77 absorption and pair production in the blob. In the simulations, the exact, 
angle-averaged synchrotron spectrum of an isotropic electron population as given by Crusius & Schlickeiser 



( |1986D is used. 

Fig. ^ shows two examples of fits to low and high flux levels of Mkn 501. The respective fit parameters 
are given in the figure captions. Inverse-Compton cooling on accretion-disk photons may be important 
close to the base of the jet even if the resulting photon spectra do not contribute significantly to the 
time-averaged emission and is thus still included in our simulations. 

Due to the non-linear nature of the model system, the fit results cannot be described in a simple, linear 
way as a function of the model parameters. We therefore construct a three-dimensional mesh of simulations 
in parameter space, with the electron density rie, the high-energy cutoff 72 and the spectral index p of the 
injected electron spectrum as parameters which are free to vary on the grid points. 

We calculate our parameter grid varying p in steps of 0.025 between 1.6 and 2.8, the total electron 
number density rie in steps of 5 cm~'^ between 10 and 120 cm^'^, and 72 for values of 2 • 10^, 3 • 10®, 5 • 10®, 
7.5 • 10®, 10^, 1.5 • 10^, 2 • 10^, 2.5 • 10^ and 3 • 10^. We constrain the range of 72 values to 72 < 3 • 10^ because 
of the kinematic limit and because around this energy, electron cooling due to triplet pair production on the 
highest-energy synchrotron photons, which is ignored in our simulations, becomes dominant over Compton 



scattering (Mastichiadis et al. 1994 , Anguelov et al. 1999 ). We find that our simulated spectra are only 



very weakly dependent on the actual value of 72. A change of 72 by a factor 3 typically results in an 
increase of the reduced of only 0.3 such that the above-mentioned restriction of 72 values has only minor 
impact on our results. In fact 72 can be regarded as constant and of the order of 10^. We point out that 
the instantaneous synchrotron spectra of individual blobs at the time of injection, which might correspond 
to short-term X-ray and TeV fiares, have their synchrotron peak at Vsy,inst. ~ 2.8 • 10® {B/G)D^2 Hz if 
p < 3, in agreement with the shift of the synchrotron peak into the hard X-ray regime during extreme 



flaring activity (e. g.. Plan et al. 1998) 



All other parameters (in particular, the magnetic fleld at the particle injection site, Bq = 0.05 G, the 



- 11 - 



low-energy cutoff of the electron spectrum 71 = 300, and the Dopplcr factor, D = 30) are fixed to values 
allowing good fits to the observed weekly averaged SEDs using our simulation code. An estimate for the 
required parameters can be found on the basis of the location of the synchrotron and SSC peaks of the 
observed broadband spectra as described below. 

Although we are assuming the injection of a single power-law distribution of ultrarelativistic electrons 
into the jet, the time-averaged radiation spectrum will be reasonably well approximated by the one produced 
by a broken power-law distribution of electrons with spectral index p below the break energy 7f,, and p + \ 
above the break energy. 7f, may be computed by setting the synchrotron cooling time scale equal to the 
dynamical time scale of jet evolution (magnetic field decay), which yields 



7b 



6.4- 10^ 



25 



^0.03 



(17) 



where = r/25, the height of the injection/acceleration site above the accretion disk is Zi = 0.03 zq. 03 pc, 
and = _Bo/(0.1 G). Thus, our model calculations will produce a time-averaged synchrotron break at 



3.4 • 10^^ S_i L>, 



30 ■ 



^ 25 



54 



Hz 



(18) 



where is an appropriate average of the magnetic field (in units of 0.1 G) over the jet evolution. For 
the purpose of these estimates, we neglect factors of (1 -I- z) ^ 1 for Mkn 501. The location of the peak of 
the SSC component will be strongly influenced by Klein-Nishina effects and will thus depend on the actual 
shape of the synchrotron spectrum, which, in turn, depends on the electron spectral index p. Considering 
these effects, Tavecchio et al. ( |1998D find 



esse ~ 76-Dg(ai,a2) 



where ai = (p — l)/2, 012 =p/2, and 

ff(ai,a2) 



exp 



ai — 1 2 (a2 — ai) 
For p = 2.5, this yields esse ^ 9-4 • 10^ D^q r25/(zo.03 ^-i)j corresponding to 



Esse 



490^^i^GeV. 



ZO.03 -S?-! 

Combining Eqs. |l^ and 21 and using the average observed e 



sy 



D: 



30 



300^g"(ai,a2) 



10 ^ and esse ^ 10^, we find 
- 0.34 



(19) 



(20) 



(21) 



(22) 



for p = 2.5 (see also Tavecchio et al. 1998). Similarly, we may use Eq. (22) of Tavecchio et al. (1998) to 
estimate 

^„2+ai^r , X i(i-Qi)/2 /2/(ai,Q;2) {vL^'^ 
BD+°'^ ^ [g{ai , aa) esse esy] W ^3 y=~ 



' sy 



SSC 



where /(ai,a2) = 1/(1 - ai) -H/ (q2 - 1 )- Using p = 2.5, {vLi,)^y ~ 6- lO'^'^ erg/s, {vL^,) 
and t^ar '^-^ 5 h (Aharonian et al. 1999a), we find 



isse 



(23) 
2-10'''' erg/s. 



B-i dH'^ > 0.34. 



(24) 
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we have D30 1 ^'^'^ ^-1 ^ 0.34. These numbers are consistent with the 
Protheroe 



1999, but are shghtly outside the allowed region of parameter 



Combining this with Eq. E 
limits found by Bednarek 

space as found by Tavecchio et al. ( 1998 ) and Kataoka et al. ( 1999 ). This is because in those papers 
the broadband spectrum is either characterized by quantities pertaining to individual outbursts or to a 
long-term quiescent state. Those parameters are not representative of the weekly averages investigated in 
this paper. 

Having constructed the three-dimensional mesh of simulations, we compare all weekly SEDs with the 
simulated spectra and find the simulation with the smallest x^. Results of this procedure are described in 
the next section. 



3. Results 

3.1. Correlation TeV-X-Ray 

The SSC model for Mkn 501 predicts a very strong correlation between the emission at the synchrotron 
peak (in soft - hard X-rays) and at the inverse-Compton peak (close to TeV 7-ray energies). We have 
derived an analytic estimate for the expected correlation, for variations of several input parameters. In 
the following discussion, unprimed quantities are measured in the co-moving frame, while a superscript * 
refers to quantities measured in the observer's frame. We assume that the time-averaged (cooling) electron 
spectrum can be described by a broken power-law, 

„ f \ / (7/7&)"^ for 71 < 7 < 7f) fr,r\ 
I (7/76) ^ for 76 < 7 < 72 

where the injection spectral index 2 < p < 3, and 7b is the break energy of the spectrum, determined by 
Eq. The normalization is given by uq w Jj^'^ / {p — 1). We are using a S approximation for the 
synchrotron spectrum: 

Lsy{e)^L,.\l^/;^\Z fo-i^^^^^ (26) 
L(e/eb) 2 for efc < e < 62, 

where e = hv/{me(?) is the dimensionless photon energy and = 2.3 • 10~^^ {B/G) 7? is the characteristic 
synchrotron energy radiated by an electron of Lorentz factor 7, . Normalizing the synchrotron luminosity to 

Lsy(xB^ J djne{-f)-/^, (27) 
71 

we have 

0.^(21)'". (28) 
P-I V76/ 

Neglecting Compton scattering events in the Klcin-Nishina regime, £7 > 3/4, we may approximate the 
SSC spectrum by 

Lsscie.)'xjde^^^n,y^^ 9(3/4-^^7^), (29) 
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where Q is the Heaviside function. The evaluation of this expressions is straightforward. Observed 
fluxes in the ASM, BATSE, and HEGRA energy ranges are calculated integrating the Doppler boosted 
synchrotron and SSC spectra, L*{e*) — L{e* /D), over the energy ranges 4 • 10^'^ < e-\sM — ^ ' 10~^, 
4 • 10-2 < e*s^^sE < 0-4, and 3 • 10^ < e^^^cflA < 6 • 10^. 

In Fig. |, we plot trajectories in the {Fbatse, Fhegra) and {Fasm i Fhegra) planes of these 
solutions, varying individual model parameter separately while all others are held constant at values 
representative of states of moderate X-ray and high-energy 7-ray fluxes. 

A variation of the electron density obviously yields a relation Fssc F^y since the synchrotron flux 
depends linearly, the SSC flux quadratically on rig. Note that this dependence may be altered due to an 
increasing 77 absorption opacity intrinsic to the source, which is not included in the analytical estimate 
(^) used to compute the HEGRA flux. A variation of the electron injection spectral index p results in a 
relation which may be approximated by Fhegra oc Fg\rpg^ and Fhegra oc F\g]yj, i. e. the dependence 
is weaker than quadratic. 

A variation of the magnetic field strength leads to more complex flux variations due to the back- 
reaction on the break Lorentz factor 7f, as a result of radiative cooling. For relatively strong magnetic fields 
{B ^ 0.3 G), the X-ray and TeV 7-ray fiuxes become anti-correlated. 

Finally, if the variability of this source were dominated by a variation of the bulk Lorentz factor F, 
the X-ray and high-energy 7-ray fiuxes would be expected to be approximately linearly correlated (the 
back-reaction on 7f, leads to a slight deviation from a strictly linear correlation), as long as the observer is 
located within the l/F beaming cone of the jet. If F increases beyond ^ l/9obs, both the X-ray and TeV 
7-ray fiuxes start to decrease with increasing F. The same quasi-linear correlation would be expected if the 
variability were due to a bending jet, i. e. a variation of 6obs- 

The empirical correlation of the TeV and the X-ray emission of Mkn 501 in 1997 has already been 
studied extensively using the CT1/CT2 and the CTS data from HEGRA and the soft X-ray data from 
RXTE ASM: Aharonian et al. ( |1999b ) find the correlation coefficient for the daily averages to be 



0.61 ± 0.057. This maximum correlation is found for zero time-lag. We examine the correlation between the 
weekly RXTE, BATSE and HEGRA points from table ||. Fig. || shows the observed correlation between the 
HEGRA and BATSE measurements, fitted with a second-order polynomial as well as with a power-law with 
index 1.4, which is the theoretical prediction if the variability is caused solely by variations of the electron 
injection spectral index p. Both fits give acceptable values for the reduced (1-1 and 1.65 respectively). 
Still, due to the large error bars, we can not confidently distinguish between these and similar correlations 
on the basis of the currently available data. In any case, there is no indication for a super-quadratic 
dependence between the X-ray and TeV fluxes, which would be inconsistent with a pure SSC model, unless 
there is a persistent quiescent level of emission above which the observed flaring behaviour is superimposed. 

Figure 1^ shows the correlation between the weekly TeV and Soft-X-ray points. The linear correlation 



coefficient is 0.59, nearly the same as found by Aharonian et al. (1999b) who compare the same data on a 
daily basis. The constant term of the linear fit is still consistent with zero. However, the reduced of 4.7 is 
too large for a good fit. This is also the case for a fitted power-law with index 1.6 which gives — 4-9. The 
systematic differences in time-coverage which are not taken into account in the determination of the error 
bars may be responsible for this. Still, the large linear correlation coefficient suggests that the correlation is 
nearly linear. 

Figure shows the correlation between the weekly Hard- X-ray and Soft-X-ray points. In this case. 
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the time coverage is the same for both instruments, only the duty cycles are different. The correlation 
coefficient is 0.53 corresponding to a 0.5 % chance probability for a linear correlation. The constant term 
of the linear fit is very well consistent with zero. This figure also illustrates the difference in the dynamical 
ranges of the variability in the soft and the hard X-ray band. At BATSE energies, which are believed to 
be near the high-energy end of the synchrotron spectrum, the variability amplitude is about 50 % larger 
than at RXTE energies. This fits nicely into the scheme that the strongest variability takes place at the 
high-energy ends of both spectral components. 



3.2. Model fit results 

Each weekly SED is compared to a three-dimensional mesh of 48 x 22 x 9 = 9504 simulations (48 
different values of p, 22 values of Tig, and 9 values of 72), selecting the simulated broadband spectrum with 
the smallest . The resulting best-fit parameters are listed in table H, along with the resulting divided 
by the number of data points. Since we do not have a continuous sequence in parameter space, the actual 
number of degrees of freedom is questionable so that we use the above quantity to assess the quality of the 
fit. 

Our best-fit parameters for each individual weekly averaged broadband spectrum are listed in Tab. ^. 
With few exceptions (MJD 50549.0, 50576.9, 50618.8, 50696.1), aU fits resulted in acceptable values. 

We find a correlation between the electron injection spectral index p and the X-ray and high-energy 
7-ray fluxes, while there is no obvious correlation with the total electron density and/or the high-energy 
cut-off 72 of the injected electron spectrum. There also appears to be a weak anti-correlation between the 
jet filling factor / c>c (Airep)"^ and the X-ray and HEGRA fluxes. This could indicate that during states of 
relatively low activity, the fluxes are dominated by a quasi-steady component from a continuous jet, while 
in high-activity state the emission is dominated by more isolated, eruptive events. However, this latter 
correlation is much less pronounced than the correlation with the electron spectral index and will need to 
be tested on the basis of future, more sensitive observations. 

In Fig. ^ the temporal variation of the best-fit values of p and Ai^ep are compared to the variations 
of the soft and hard X-ray fiuxes and the 1.5 TeV flux. The correlation between the hard X-ray and TeV 
7-ray fluxes with the injection spectral index is illustrated in figure |l^. 

Furthermore, we show in figure |ll| the correlation between p and the positions of the peaks in the 
synchrotron and the inverse Compton component of the SED. The peak positions are not fit parameters 
but were determined by finding the local maxima in the weekly SEDs. We find that, in our model, there is 
a strong correlation between the peak positions and p such that these parameters can be regarded as nearly 
identical. However, while the peak positions are directly observable, the electron spectral index is the more 
fundamental quantity. 

Our results indicate that medium-timescale high activity states in X-rays and high-energy 7-rays are 
consistent with a hardening of the electron spectrum injected at the base of the jet. As pointed out in the 
previous subsection, a pure SSC model in which the TeV 7-ray fiux is strongly infiuenced by Klein-Nishina 
effects, predicts that the HEGRA and both the soft and hard X-ray fluxes should roughly be correlated 
by power-laws Fhegra oc with 1.4 <J (5 <^ 1.6 in high flux-level states, in which the contribution of a 
possible quasi-stationary radiation component is small. The data available for this study are consistent with 
this but do not allow a clear distinction between different variability mechanisms, and future observations 
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with increased sensitivity, in particular at multi-GeV to TeV energies, are needed in order to test this 
prediction. 



4. Summary and Conclusions 

We have presented broadband spectra of the extreme HBL Mkn 501 during its high state in 1997, 
including radio, optical, soft and hard X-ray, and TeV 7-ray observations. In this study we concentrated 
on the medium-timescale variability, using weekly averaged SEDs. We confirmed the strong correlation 
between the TeV 7-ray flux and the hard X-ray flux. This correlation was found to be non-linear and could 
be fitted with a second-order polynomial, in agreement with the expectation of an SSC dominated leptonic 
jet model, if the flux variations are related to fluctuations of the electron density in the jet and/or the 
spectral index of the electron spectrum at the time of injection into the jet. 

The weekly averaged SEDs were fitted with a leptonic jet model, strongly dominated by the SSC 
process. With a few exceptions, this model yielded acceptable fits to the observed broadband spectra. The 
observed spectral variability of Mkn 501 could be explained mainly by variations of the electron spectral 
index. No clear correlation between the maximum electron energy and the hard X-ray and TeV 7-ray fluxes 
on the 1-week timescale was found, in contrast to the short-term variability of Mkn 501. Plan et al. ( |1998| ) 
have shown that the intraday variability of this object is most probably related to an increase of 72, leading 
to pronounced flares in hard X-rays, most probably on the synchrotron cooling timescale which is most 
likely a few hours and thus much shorter than the 1-week timescale considered in this paper. Our result 
indicates that such synchrotron flares are isolated events and are at most weakly correlated to the activity 
of the source on the 1-week timescale. 

Our result that the flaring behaviour on intermediate timescales is consistent with a hardening of the 
electron spectrum is in contrast to the flaring characteristics observed in quasars. Recently, Mukherjee et 
al. ( 1999| ) have investigated all available broadband data on the very luminous flat-spectrum radio quasar 



(FSRQ) PKS 0528-1-134, and found that its flaring behaviour is consistent with an increasing contribution of 
the external inverse- Compton component during flares, possibly related to an increase in the bulk Lorentz 
factor. The flts to the SEDs of PKS 0528-f 134 required that the average energy of relativistic electrons in 
the jet shifts towards lower values during flares, in contrast to the results found for Mkn 501. As pointed 



out by Bottcher (1999), this implies that the synchrotron peak is expected to shift towards lower frequencies 
during flares of FSRQs, while Mkn 501 and Mkn 421 show clear evidence for a shift of the synchrotron peak 
to higher frequencies. 

We point out that in the present study the magnetic fleld along the jet and the bulk Lorentz factor 
of individual blobs were fixed, so that we cannot confidently rule out variations of the Doppler factor, 
accompanied by appropriate changes of the electron injection spectrum, as the flaring mechanism for Mkn 
501. However, the very moderate variability at optical frequencies, as observed in Mkn 501, leads us to 
consider this flaring mechanism less likely in this object since it would require a peculiar conspiracy between 
the Doppler factor and the electron spectrum to keep the optical flux at an approximately constant level. 
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Table 1: The instruments which contributed data to this paper. 



Instrument 


energy range 


comment 




(Hz) 


(cV) 




Metsahovi Radio Telescope 


22 X 10® 


9 X 10"^ 


A = 1.4 cm 


Nordic Optical Telescope 


4.4 X 10" - 6.2 X 10" 


1.8 - 2.6 


2.5 m mirror, filters: BVRI 


Tuorla Observatory 




J5 


1.0 m mirror, filters: V 


Osscrvatorio di Torino 




): 


1.0 m mirror, filters: BVR 


Osservatorio di Perugia 




J5 


0.4 m mirror, filters: VRI 


RXTE ASM 


4.8 X 10^^ - 2.4 X 10^* 


2 X 10^ - 1 X 10* 


(see e.g. Levine et al. 1996) 


BATSE 


4.8 X 10^* - 4.8 X 10^® 


2 X 10* - 2 X 10^ 


occultation measurement 


HEGRA CTl 


3.6 X 10^® - w 7.3 X 10^'' 


1.5 X 10^^ - 3 X 10" 


Aharonian et al. (1999b) 


HEGRA CT System 


1.9 X 10^® - w 1.2 X 10^* 


8 X 10" - w 5 X 10" 


Aharonian et al. (1999a) 



Table 2: The weekly i/F^ datapoints from observations of Mkn 501 used for the model fits presented in this 
paper. All entries are in units of 10~^^ erg cm~^ s~^. Entries exactly equal to are those where no data 

are available. This is mainly radio data. Upper limits arc indicated by a "<" symbol and arc calculated 
with 90% confidence kn'cl. TIk^ fluxc^s at 0.8 TcV and 10 TcV iu~c cxtrapolatcxi from tlic^ flux mcasiu'cmcnt 
at 1.5 ToV. 



MJD 


Radio 


Optical 




Soft X-ray 


Hard X-ray 


(0.8 TcV) 


1.5 TcV 




(10 TcV) 


50517.199 


0.02244 ± 0.00104 


7.068 ± 0, 


.124 


11.1 ±3.1 


26.73 ± 8.5 


8.7 ± 1.9 


8.2 ± 1, 


.9 


2, 


,41 ±0.65 


50521.906 


0.02387 ± 0.00156 


7.148 ±0, 


.069 


10.0 ±3.1 


1.67 ±9.5 


10.6 ± 1.5 


9.9 ± 1 


.5 


2, 


,94 ±0.61 


50526.34 


0.02200 ± 0.0011 


7.275 ± 0, 


,047 


15.8 ±3.3 


14.33 ±6.2 


7.4 ± 1.4 


6.9 ± 1 


.3 


2, 


,05 ± 0.49 


50536.715 


0.02420 ± 0.0011 


7.394 ± 0, 


.120 


17.4 ±3.5 


8.38 ± 8.4 


5.6 ± 1.9 


5.3 ± 1, 


.9 


1, 


.57 ±0.59 


50541.262 


0.02237 ± 0.00064 


7.607 ± 0, 


.077 


11.9 ±3.0 


< 8.4 


7.8 ±1.7 


7.3 ±1 


.7 


2, 


.17 ±0.58 


50549.012 





7.769 ± 0, 


.063 


15.3 ±3.8 


45.7 ±7.4 


13.4 ± 1.4 


12.6 ± 1 


.5 


3, 


,72 ± 0.69 


50556.32 


0.02508 ± 0.00198 


7.708 ± 0, 


.109 


20.9 ± 3.9 


36.38 ± 6.9 


17.6 ± 2.4 


16.5 ±2 


.4 


4, 


.88 ± 0.99 


50564.566 


0.02382 ± 0.00052 


7.948 ± 0, 


.093 


21.9 ±3.8 


16.79 ±8.6 


7.9 ±1.2 


7.4 ±1 


.2 


2, 


.18 ±0.47 


50570.254 


0.02420 ± 0.00095 


7.941 ± 0, 


,075 


20.1 ±4.3 


16.34 ± 5.8 


6.8 ±0.8 


6.4 ±0 


.8 


1, 


,90 ± 0.36 


50576.934 


0.02409 ± 0.00104 


8.071 ± 0, 


.188 


20.6 ± 3.9 


49.87 ± 8.2 


16.2 ± 1.5 


15.2 ± 1 


.7 


4, 


.49 ± 0.80 


50583.305 





7.921 ± 0, 


.317 


18.1 ±3.2 


37.24 ± 6.6 


16.9 ±1.9 


15.8 ±2 


.0 


4, 


.68 ± 0.89 


50592.621 





7.829 ± 0, 


,153 


16.9 ±3.2 


13.03 ± 7.7 


5.7 ± 1.8 


5.3 ±1 


.7 


1, 


,58 ±0.55 


50600.797 





7.393 ± 0, 


.188 


28.3 ±4.7 


29.0 ± 10.3 


6.9 ± 1.3 


6.5 ± 1, 


.2 


1, 


,92 ± 0.46 


50604.852 


0.02574 ± 0.00132 







35.7 ±5.6 


43.89 ± 9.0 


17.2 ± 1.8 


16.2 ± 1 


.9 


4, 


,78 ±0.88 


50611.523 





7.029 ± 0. 


.317 


26.0 ±5.4 


37.29 ± 11 


13.3 ± 1.5 


12.5 ± 1 


.6 


3, 


,68 ±0.69 


50618.848 










35.5 ± 5.4 


52.0 ± 7.4 


3.0 ±2.0 


2.8 ± 1, 


.8 


0, 


,83 ±0.56 


50626.785 


0.02420 ± 0.00154 







37.2 ± 6.2 


39.19 ±6.1 


20.7 ±2.1 


20.0 ± 2 


.3 


5, 


,90 ± 1.07 


50634.707 


0.02189 ± 0.00148 


7.327 ± 0, 


,070 


16.1 ±3.7 


10.61 ± 7.4 


3.3 ±0.6 


3.1 ±0 


.5 


0, 


,92 ±0.21 


50640.742 





7.271 ± 0, 


.138 


28.8 ±4.7 


36.5 ± 8.6 


18.2 ± 1.8 


17.1 ±2 


.0 


5, 


,06 ±0.92 


50650.969 





7.859 ± 0, 


.142 


26.0 ±4.5 


45.82 ± 7.9 


9.4 ± 1.8 


8.8 ± 1, 


.7 


2, 


,60 ± 0.63 


50654.59 










19.4 ±3.7 


28.48 ± 9.7 


8.5 ± 1.3 


8.0 ± 1, 


.3 


2, 


,36 ±0.52 


50661.887 


0.01672 ± 0.00132 


7.226 ±0, 


.317 


10.4 ±4.6 


24.66 ± 7.7 


9.2 ± 1.0 


8.7 ± 1, 


.1 


2, 


,57 ±0.48 


50669.164 


0.02706 ± 0.00154 


7.567 ± 0, 


.217 


14.0 ±4.8 


36.78 ± 15 


11.3 ± 1.5 


10.6 ± 1 


.5 


3, 


,12 ±0.63 


50672.938 





7.244 ± 0, 


,160 


20.4 ±4.7 


30.82 ±6.8 


4.7 ± 1.9 


4.5 ± 1 


.8 


1, 


,32 ±0.57 


50684.598 





7.484 ± 0, 


.177 


24.2 ± 5.8 


38.19 ±5.8 


14.6 ± 1.7 


13.7 ± 1 


.8 


4, 


,04 ± 0.77 


50689.281 





7.162 ±0. 


.094 


16.6 ± 3.4 


12.31 ± 7.9 


6.6 ± 1.1 


6.2 ± 1 


.2 


1, 


,84 ± 0.42 


50696.059 










13.5 ±3.3 


45.06 ±8.9 


7.1 ± 1.1 


6.7 ± 1 


.1 


1, 


,98 ± 0.43 


50703.445 










20.6 ± 3.7 


13.63 ± 8.8 


9.4 ±2.8 


8.8 ±2, 


.7 


2, 


,61 ±0.87 
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Table 3: Best-fit parameters of a pure SSC jet model to the weekly SEDs, and the quality of the fit, 
indicated by / (no. of data points). The is calculated without taking into account the radio point. 
Periods marked with a have less than 3 HEGRA points. Vg and Vic arc respectively the positions of the 
peak in the synchrotron and inverse compton component of the SED. They were determined by finding the 
two local maxima in each fitted SEP. 



MJD-50000 72/10'^ 




P 


rie 


At [10* s] 


log(i^s/Hz) 


log(i/. 


c/Hz) 


1 2 
red. X 


no. data points 


517.199 


3.0 


2, 


.425 


30 


3.443 


17.23 


25, 


.04 


0.821 


7 


521.906 


3.0 


2, 


.450 


65 


6.189 


16.34 


24, 


.74 


0.395 


7 


526.340 


3.0 


2 


.425 


20 


2.308 


17.51 


25, 


.26 


0.291 


7 


536.715 * 


3.0 


2, 


.450 


15 


1.605 


17.54 


25, 


.11 


0.254 


7 


541.262 


3.0 


2, 


.500 


50 


4.071 


16.51 


24, 


.77 


0.539 


7 


549.012 


3.0 


2, 


.350 


25 


3.285 


18.00 


25, 


.30 


3.264 


6 


556.320 


3.0 


2, 


.300 


25 


3.798 


18.53 


25, 


.40 


0.806 


7 


564.566 


3.0 


2, 


.400 


15 


1.737 


17.76 


25, 


.30 


0.326 


7 


570.254 


3.0 


2, 


.425 


15 


1.614 


17.62 


25, 


.28 


0.495 


7 


576.934 


3.0 


2, 


.300 


20 


2.991 


18.53 


25, 


.49 


2.039 


7 


583.305 


3.0 


2, 


.325 


30 


4.048 


18.18 


25, 


.32 


1.760 


6 


592.621 * 


3.0 


2, 


.475 


20 


1.849 


17.34 


25, 


.08 


0.241 


6 


600.797 


3.0 


2, 


.375 


10 


1.363 


18.00 


25, 


.49 


0.412 


6 


604.582 


2.5 


2, 


.200 


10 


2.476 


19.00 


25, 


.74 


0.309 


6 


611.523 


3.0 


2, 


.300 


15 


2.609 


18.42 


25, 


.52 


0.483 


6 


618.848 * 


2.0 


2, 


.400 


10 


1.146 


17.81 


25, 


.40 


5.177 


5 


626.785 


3.0 


2, 


.225 


15 


2.989 


19.04 


25, 


.58 


0.140 


6 


634.707 


3.0 


2, 


.500 


10 


0.940 


17.51 


25, 


.08 


0.522 


7 


640.742 


3.0 


2, 


.250 


15 


2.941 


18.90 


25, 


.54 


0.354 


6 


650.969 * 


3.0 


2, 


.325 


10 


1.500 


18.23 


25 


.53 


1.326 


6 


654.590 


3.0 


2 


.375 


15 


2.039 


17.78 


25, 


.40 


0.560 


5 


661.887 


3.0 


2, 


.425 


30 


3.361 


17.26 


25, 


.15 


0.959 


7 


669.164 


3.0 


2, 


.375 


25 


3.151 


17.72 


25, 


.28 


0.698 


7 


672.938 * 


3.0 


2, 


.400 


10 


1.291 


17.78 


25, 


.38 


0.706 


6 


684.598 


3.0 


2, 


.275 


15 


2.676 


18.64 


25, 


.53 


0.849 


6 


689.281 


3.0 


2, 


.450 


20 


2.175 


17.36 


25 


.08 


0.378 


7 


696.059 


3.0 


2 


.425 


20 


2.319 


17.51 


25 


.15 


2.996 


5 


703.445 


3.0 


2, 


.375 


20 


2.629 


17.76 


25, 


.30 


0.197 


5 


Fit parameters: 



72 = electron spectrum high-energy cutoff 

p = electron spectral index, 71(7) oc 7"'' 

Ue [cm~^] = electron density 

At [s] = blob ejection events repetition time scale (normalization) 
Fixed parameters: 

Zi = 0.03 pc (injection height of blob) 

Mbh = 10® Mo (mass of central black hole) 

Ld = 5 X 10** ergs" ^ (isotropic accretion disk luminosity) 

Rg = 3 X lO'^^ cm (blob radius in the comoving frame) 

71 = 200 (low-energy cutoff of electron sp.) 

r = 25 (bulk Lorentz factor) 

D = 30 (Doppler factor) 

Stmin = 3336 s (contracted blob crossing time) 

B = 0.05 G (magnetic field) 

Ho = 75 kms"^Mpc"^ (Hubble constant) 

go = 0.5 (deceleration parameter) 

z = 0.034 (cosmological redshift) 

= 133 Mpc (luminosity distance) 
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Fig. 1. — The lightcurve data which were used in this paper to construct weekly spectral energy distributions. 
See section for references. 
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Fig. 2. — The correction factor k for the second RXTE ASM energy bin derived from the pubhcly available 
Crab Nebula data taken by the detector between MJD 50510 and 50710. See text. 
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Fig. 3. — The index a of the differential photon spectrum in the energy range 3.0 - 12.1 keV derived from 
the publicly available RXTE ASM data for Mkn 501 and the Crab Nebula in weekly time bins. The lines 
represent fits of constant functions. For Mkn 501, the reduced of the fit is 0.93, for the Crab it is 0.95 . 
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v[Hz] 

Fig. 4. — Model fits to the weekly averaged broadband spectra of Mkn 501 for the periods centered on 
MJD 50564.566 (low flux state; filled triangles and dot-dashed curve) and MJD 50626.785 (high flux state; 
filled circles and solid curve). Model parameters for MJD50564.566: 71 = 500, 72 = 3 • lO'^, p = 2.400, 
He = 15 cm-3, Atrep = 1-74 • 10^ s, B = 0.05 G, r = 25, i?^ = 3 • 10^^ cm, D = 30. Model parameters for 
MJD 50626.785: Same as for the low state, except p = 2.225, Atrep = 3.0 • 10^ s. 



- 25 - 



10 



10 



w 
3 



10° 



< 

^2 7 

^ 10 



10" 









/ '/ : 


var. oT p 






/' / 
/ / 


Var. of B 






A 










r Var. of 


"e = 




y/ p = 2.2 : 






















r = 50 






8 = 0.1 G 




// 






_> 






r = 20 




\ 








, B = 1 G 






/v' 


p = 2.5 


B = 0.01 G 


' ' ' ' ■•■'/^ 














■/••'-""' 






50 cm ^ 


r=io 










/ \ 






.>-"""' 


/ /' ^P = 


2.7 






/ 








n. = 1 cm" 

S 












10" 




10-^ 


10° 



Fasm (arbitr. units) 



10 



10 



3 



CO 



lO'' 



10° 



< 

DC 

^ 7 

^ 10 



10" 







/ 


■ var. of p 




/' 

/' 


var. of B 




/' 


Var. of r 




/ /y 


r Var. of 




/' yA 




= 200 cm'' 


/' yA / 






















B = 0.1 G 






r = 20 ^~>\ 


" r=50 






, ^B = 1 G 












^p = 2.5 


B = 0.01 G 










50 cm"^ 










/X^P = 2.7 




V=10 


/ \ 






/' ^n^ = 20cm"^ 










10-^ 


10-^ 

Fbatse (arbitr. units) 


10° 



Fig. 5. — Correlation between the ASM and HEGRA fluxes (a) and between the BATSE and HEGRA 
fliixes (b) according to our analytical approximation. Standard model parameters arc 71 = 300, B = 0.1 G, 
rie = 100 cm~^, p = 2.5, 72 = 10^. For each curve, one parameter is varied, while the others are fixed to the 
above values. The curves are labelled by a few representative values of the varying parameter. 
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BATSE i/F, (10 "^^ erg cm'' s"') 

Fig. 6. — The correlation between HEGRA TeV and the BATSE X-ray flux from Mkn 501. In order to avoid 
effects from poor time-coverage, points obtained from less than three independent HEGRA measurements 
were excluded. The excluded points are shown as open circles. The fit of a second-order polynomial to the 
remaining points (filled circles) yields y = (5.4 ± 0.6) + (-0.29 ± O.G8)a; + (0.017 ± 0.0035)a;2 (solid line), 
where y = HEGRA flux at 1.5 TeV in units of 10"" erg cm^^ g-i^ ^nd x = BATSE flux at 36.4 keV in 
the same units, and results in a reduced of 1.14. With the open circle points included the reduced 
increases to 1.54. The dashed line is a linear fit to the filled circle points: y = (—0.3 ± 1.7) + (0.40 ± 0.06)a;, 
reduced = 1-42 (2.6 with the open circle points). The dot-dashed hue is a fit of the function y = ax^-'^ -\- b. 
It results in a = 0.063 ± 0.0084, 6 = 4.2 ± 0.72, and a reduced of 1.65 (2.7 with the open circle points). 
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Fig. 7. The correlation between the HEGRA TeV and the RXTE ASM X-ray flux from Mkn 501. In 
order to reduce effects from poor time-coverage, points obtained from less than three independent HEGRA 
measurements were excluded. The excluded points are shown as open circles. The fit of a linear function 
to the remaining points (filled circles) yields y = (-2.6 ± 1.40) + (0.62 ± 0.07)a;, where y = HEGRA flux at 
1.5 TeV in units of 10^^^ erg cm^^ s~^, and x = RXTE flux at 5.2 keV in the same units, and results in a 
reduced of 4.7. Only statistical errors were taken into account. The linear correlation coefficient is 0.59. 
The dot-dashed line is a fit of the function y = ax'^-^ + b which results in a = 0.062 ± 0.0075, b = 2.28 ± 0.80, 
and a reduced of 4.9 . 
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Fig. 8.— The correlation between the BATSE hard X-ray and the RXTE ASM soft X-ray flux from Mkn 501. 
Here the time coverage is not systematically different, only the BATSE duty cycle is much lower than that 
of the RXTE ASM. So all points are used in the linear fit. The points where HEGRA had bad time coverage 
are are still marked as open circles for comparison with the other figures. The fit of a linear function yields 
y = (-0.1 ± 5.3) + (1.36 ± 0.2A)x, where y = BATSE flux at 36.4 keV in units of lO'^^ erg cm'^ s-\ and 
X = RXTE flux at 5.2 keV in the same units, and results in a reduced of 2.3. Only statistical errors were 
taken into account. The linear correlation coefficient is 0.53. 
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Fig. 9. Best-fit values of tlie blob ejection repetition time At^ep and the electron injection spectral index 
p for the pure SSC model compared to the RXTE ASM, BATSE and HEGRA 1.5 TeV hght curves. 
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Fig. 10. — The correlation between the fit parameter spectral index, p, and the HEGRA TeV (filled 
circles) and BATSE hard X-ray (stars) flux from Mkn 501. The fit of a linear function yields for HEGRA 
j5= (2.57 ±0.018) + (-0.019 ±0.0018)a;, where a; = HEGRA fiux at 1.5 TeV in 10"" erg cm'^ 8"^ reduced 
= 1.18: the linear correlation coefficient is -0.89. For BATSE: p = (2.56 ± 0.034) + (-0.0072 ± 0.0011).t, 
where x = BATSE flux at 36.4 keV in 10^^^ erg cm~^ s~^; reduced = 1.2; the linear correlation coefficient 
is -0.76. 
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Fig. 11. — The correlation between the fit parameter spectral index, p, and the positions of the peak in 

the synchrotron (bottom plot) and inverse Compton (top plot) component of the weekly SEDs. The fit of 
a linear fimction yields log(i^s/Hz) = (39.11 ± 0.68) + (— 8.97± 0.29)p (linear correlation coefficient —0.93) 
and log(z/ic/Hz) = (31.84 ± 0.38) + (-2.77 ± 0.16)p (linear correlation coefhcient -0.88). 



